Optimum pulse shapes for stimulated Raman adiabatic passage 
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Stimulated Raman adiabatic passage (STIRAP), driven with pulses of optimum shape and delay 
has the potential of reaching fidelities high enough to make it suitable for fault-tolerant quantum 
information processing. The optimum pulse shapes are obtained upon reduction of STIRAP to 
effective two-state systems. We use the Dykhne-Davis-Pechukas (DDP) method to minimize nona- 
diabatic transitions and to maximize the fidelity of STIRAP. This results in a particular relation 
between the pulse shapes of the two fields driving the Raman process. The DDP-optimized version 
of STIRAP maintains its robustness against variations in the pulse intensities and durations, the 
single-photon detuning and possible losses from the intermediate state. 

PACS numbers: 03.65.Ge, 32.80.Bx, 34.70.+e, 42.50.Vk 



I. INTRODUCTION 

Stimulated Raman adiabatic passage (STIRAP) is a 
well established and widely used technique for coherent 
population transfer in atoms and molecules [l[. STI- 
RAP uses two delayed but partially overlapping laser 
pulses, pump and Stokes, which drive a three-state A- 
system ipi — > ■02 ^ ^3- The STIRAP technique trans- 
fers the population adiabatically from the initially pop- 
ulated state tjji to the target state ip^. If the pulses are 
ordered counterintuitively, i.e. the Stokes pulse precedes 
the pump pulse, two-photon resonance is maintained, and 
adiabatic evolution is enforced, then complete population 
transfer from tpi to ip3 occurs. Throughout this process, 
no population is placed in the (possibly lossy) interme- 
diate state Various aspects of STIRAP have been 
subjects of intense research, both theoretically and ex- 
perimentally 0. 

Because STIRAP is an adiabatic technique it is insen- 
sitive to small to moderate variations in most of the ex- 
perimental parameters, such as pulse amplitudes, widths, 
delay, and single-photon detuning. A particularly re- 
markable and very useful feature of STIRAP is its in- 
sensitivity to the properties of the intermediate state "02 • 
For instance, STIRAP has been demonstrated with high 
efficiency even for interaction durations exceeding the 
lifetime of ip2 by a factor of 100 [2,]. For these reasons 
STIRAP is a very attractive technique for quantum in- 
formation processing (QIP) [1, However, it is widely 
recognized that QIP requires very high fidelities, with 
the admissible error of gate operations being below 10~^ 
for a reliable quantum processor 0] . Such an extreme 
accuracy has not yet been demonstrated for STIRAP, 
as an accuracy of 90-95% was sufficient for most tradi- 
tional applications. When trying to increase the fidelity 
beyond this number one faces various obstacles related 
mainly to nonadiabatic transitions. Being an adiabatic 
technique, STIRAP reaches an efficiency of unity only 
in the adiabatic limit; however, the latter is approached 
only asymptotically as the pulse areas increase. For QIP, 



the pulse areas needed are so large that they may violate 
various restrictions of a real experiment. 



In this paper we propose how to achieve an ultrahigh 
fidelity in STIRAP, and thus make it fully suitable for 
QIP by suitably shaped pulses. We utilize a recent idea 
of Guerin et al. Q who applied the well-known Dykhne- 
Davis-Pechukas (DDP) method Q to optimize adiabatic 
passage in a two-state system. In order to adapt this ap- 
proach to STIRAP, we reduce the three-level Raman sys- 
tem to effective two-state systems in two limits: on exact 
resonance and for large single-photon detuning. The op- 
timization, which minimizes the nonadiabatic transitions 
and maximizes the fidelity, leads to a particular relation 
between the pulse shapes of the driving pump and Stokes 
fields. 



It should be noted that a fidelity of unity can also be 
achieved by a resonant 7r-pulse in a two-state transition. 
However, resonant techniques suffer from their sensitiv- 
ity to parameter fluctuations. The optimized version of 
STIRAP presented here features both a very high fidelity 
and a robustness against variations in the intensities and 
the single-photon detuning. 



This paper is organized as follows: In section [IT] we re- 
view the DDP method and the optimization of two-state 
adiabatic passage. Then we extend this idea to STIRAP 
in section IIIII and discuss examples in section IIVI In 
section |V] we extend these ideas to fractional STIRAP 
(f-STIRAP), which creates a coherent superposition of 
Tpi and ■03. We summarize the results in the concluding 
section. 
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II. OPTIMIZATION OF ADIABATIC PASSAGE 
BETWEEN TWO STATES 

A. Dykhne-Davis-Pechukas (DDP) approximation 

The probability amplitudes in a two-state system 
a{t) — [ai{t), a2{t)\^ satisfy the Schrodinger equation, 



(1) 



where the Hamiltonian in the rotating-wave approxima- 
tion (RWA) reads @ 



n{t) 
n{t) 2A{t) 



(2) 



The detuning A = luq — lo is the difference between the 
transition frequency ujq and the carrier laser frequency 
uj. The time- varying Rabi frequency = \dE{t)\/h 
describes the laser-atom interaction, where d is the elec- 
tric dipole moment for the ipi ^ ip2 transition and E{t) 
is the laser electric field envelope. 

A very accurate technique for deriving the transition 
probability in the near-adiabatic regime is the Dykhne- 
Davis-Pechukas (DDP) approximation [8]. The DDP for- 
mula gives the following expression for the probability for 
nonadiabatic transitions 



P 



where 



-2ImD(to) 



to 



D{ta) = / e{t)dt 



(3) 



(4) 



is an integral over the splitting e{t) — y/n{ty + A{ty 
of the eigenenergies of the Hamiltonian The point 
to (the transition point) is defined as the (complex) zero 
of the quasienergy splitting, e{to) = 0, which lies in the 
upper half of the complex f-plane (i.e., with Imto > 0). 
Equation ^ gives the correct asymptotic probability for 
nonadiabatic transitions provided: (i) the quasienergy 
splitting e{t) does not vanish for real t, including ±oo; 
(ii) e{t) is analytic and single-valued at least through- 
out a region of the complex t-plane that includes the 
region from the real axis to the transition point to; (iii) 
the transition point to is well separated from the other 
quasienergy zero points (if any), and from possible sin- 
gularities; (iv) there exists a level (or Stokes) line defined 
by lmD{t) = ImZ)(to), which extends from —oo to +oo 
and passes through to- 

For the case of multiple zero points in the upper t- 
plane, Eq. ([3]) can be generalized to include the contri- 
butions from all these N zero points tk as 



P ; 



(5) 



where = 4i lim {t—tk)'d{t) \ usually F^ = 1 or —1. Here 



'd{t) accounts for the nonadiabatic coupling between the 
adiabatic states, with ^{t) 



1 tan-il7(t)/A(t). 



B. Optimization of two-state adiabatic passage 

Guerin et al. [7] have used the DDP method to opti- 
mize the adiabatic passage between two states in a very 
simple manner. Assuming that the probability for nona- 
diabatic losses is solely due to the transition points t^, 
they have proposed to suppress these altogether by choos- 
ing the Rabi frequency ri(t) and the detuning A(t) such 
that there are no transition points. This condition is ob- 
viously fulfilled if the quasienergy splitting is constant, 



£{t) = ^Jn{tY + A{tY = const. 



(6) 



For example, this condition is fulfilled for a detuning and 
Rabi frequency defined as 

A{t) = So cos [nfit)] , n{t) = eo sin [^/(t)] , (7a) 
- /(-oo) ^ fit) ^ /(oo) = 1, (7b) 

with f{t) being an arbitrary monotonically increasing 
function with the above property. Condition ([6]) is not 
the only possible condition for adiabatic optimization, 
but it is the simplest one 



III. OPTIMIZATION OF STIRAP 



STIRAP 



The probability amplitudes of the three states in STI- 



RAP c{t) 
equation. 



[ci{t),C2(t),C3{t)] satisfy the Schrodinger 



H(t)c(t), 



(8) 



where the STIRAP Hamiltonian within the rotating- 
wave approximation (RWA) reads [S] 



H(t) = i 



flpit) 

np{t) 2 A iisit) 
risit) 



(9) 



The time- varying Rabi frequencies flp(t) and ^s{t) de- 
scribe the couplings between the intermediate state V'2 
and, respectively, the initial state ipi and the target final 
state 11^3. STIRAP is easily explained with the so-called 
dark state (pd{t), which is a zero-eigenvalue eigenstate of 
H(t), 



^d[t) = Vi cos'd{t) - ipz sini?(t). 



(10) 



where the time-dependent mixing angle i?(t) is defined 
as 



i9(t) = tan- 



-1 ^pjt) 
^sit)' 



(11) 



The pulses in STIRAP are ordered counterintuitively, 
i.e., the Stokes pulse precedes the pump pulse. 



njt) 

hm = 0, ..... ^ , , 

t^-00il^{t) t^ + Q0ilp(t) 



hm 241=0. 



(12) 
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Then ^2^^* ^J(i) t—>oc ^ a,nd therefore the dark state 
^Pdify connects adiabatically states '01 and "03, 



01 



'03 



(13) 



Thus, if the evolution is adiabatic then the population 
passes from state 0i to state 03. Moreover, because the 
dark state does not contain a contribution from the (pos- 
sibly lossy) intermediate state 0)2 7 the properties of the 
latter are less important. 

When the evolution is not completely adiabatic, the 
population transfer -01 — > 0^3 might be incomplete. More- 
over, the intermediate state receives some transient pop- 
ulation, which may either be lost if state 0'2 decays to 
other states, or lead to decoherence if it decays back into 
0)1 or 0)3. For STIRAP to be a viable tool for quan- 
tum computing, nonadiabatic transitions have to be re- 
duced below the fault-tolerance limit of < 10~^ 0, [1]. 
We can estimate the required resources for STIRAP to 
reach such a fidelity as follows. The adiabatic condi- 
tion for STIRAP (for A = 0) demands large pulse ar- 
eas A^,s = J^rip^s(i)di » 1. This global condition is 
derived from the local adiabatic condition, which reads 

il(t) 3> ^9(t) ■ The probability for nonadiabatic transi- 
tions in the perturbative limit is Pna(i) ~ ^^^(i)/f^^(i); it 
measures the population that escapes to other adiabatic 
states and reduces the fidelity. The infidelity is therefore 
1 _ Pg ^ 1/Ap The fault-tolerance QIP limit therefore 
requires Ap^s ^ 100. In fact, this very rough estimate 
neglects various details^ such as the peculiarities of the 
nonadiabatic coupling "dit), and in reality the condition 
for the pulse areas is more restrictive. 

In the following, we show that an optimized version of 
STIRAP can reach the fault-tolerance QIP limit by using 
much smaller pulse areas. In order to optimize STIRAP 
we use the same ideas as in the two-state adiabatic pas- 
sage optimization by Guerin et al. [7|- To this end we 
make use of the reduction of STIRAP to effective two- 
state problems. 



B. Effective two-state systems 

1. Single-photon resonance 

On single-photon resonance (A = 0) the three-state 
system is reduced to an effective two-state system, with 
a detuning Qs{t) and a coupling Qp{t) [Tol. IT]|. 



(14) 
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'hit)' 




. b2{t) _ 
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where the probability amplitudes &i,2(i) are related to 
ci,2,3(i) as follows Ha, [H 

ci{t) = 2Re[hl{t)h2{t)]smd{t) 

+ (|6i(t)p-|62(t)ncos^(t), (15a) 
C2{t) = 2ilm[6*(t)62(t)], (15b) 
C3(t) = 2Re[6t(t)62(i)]cosi9(t) 

-(|6i(i)P-|^'2(i)nsin^(i). (15c) 

Because we have for STIRAP z9(— oo) = and -d^oo) — 
7r/2, the initial condition ci(— oo) — 1 demands the con- 
dition |6i(— oo)| = 1 in the effective two-state system. 
The final-state population in STIRAP is 



|C3( + 00)P= [|fel( + 00)|2-|62( + ^)| 



(16) 



Consequently, an optimized adiabatic evolution in the 
two-state system (|14l) implies optimized STIRAP. Hence 
applied to STIRAP, the two-state optimum condition ^ 
simply yields 



const. 



(17) 



In other words, the rms Rabi frequency should be con- 
stant. Again, as in the two-state optimization, this is 
not the only possible optimization condition but it is the 
simplest one. 



2. Far-off-resonance fields 

For large single-photon detuning A(t), the interme- 
diate state can be eliminated adiabatically by setting 
C2{t) = in Eq. ©. We thus obtain 



d 




_ 1 






" 2 



-Aeff(i) n,s{t) 

n,s{t) Aeff(t) 



Cl(t) 



(18) 



where the effective Rabi frequency flcsit) and detuning 
Aeff(<) are defined as 



f^cff(i) 
Acff(t) 



2A{t) ' 
4A(t) 



(19a) 
(19b) 



The two-state condition for optimal adiabatic passage 
now reads 



f)eff(i)' + Aeff(t)' 



npjtf + nsjty 

4A(i) 



const. (20) 



For constant A this condition is identical to the one 
we found on resonance, Eq. (jl7p . that is it requires a 
constant rms Rabi frequency fl. We point out that, due 
to the identical conditions, on and off single-photon res- 
onance, optimization is ensured over a very wide range 
of single-photon detunings. 
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IV. OPTIMIZATION OF STIRAP: EXAMPLES 
A. Pulse shapes 

Conditions (fT7|) and ((20|l suggest the following param- 
eterization of the pump and Stokes fields 



ilpit) = nsin -fit) 



Q cos 



fit) 



(21a) 
(21b) 



where f{t) is an arbitrary monotonically increasing func- 
tion, = fi-oo) S fit) S /(oo) = 1. Viewed mathe- 
matically, ftpit) and Qsit) define an adiabatic path, for 
which the nonadiabatic correction given by the DDP for- 
mula vanishes, which leads to an optimal adiabatic fol- 
lowing of the dark state. 

The exact fulfillment of condition (|17p requires a con- 
stant and hence, infinite pulse areas. This unphysical 
condition can be overcome by using a "mask" function 



npit) = noFit)sm[-fit) 
Qsit) = r!oF(t)cos [|/(i) 



(22a) 
(22b) 



Then the rms Rabi frequency becomes time- dependent, 
Hit) — riaFit). This replacement does not necessar- 
ily violate the optimization condition (|17p because DDP 
transition points may still be absent, e.g. if the mask 
function has a suitable shape, such as Gaussian. Still, 
a pulse-shaped mask Fit) violates the DDP conditions 
because the eigenenergy splitting e(i) = fioi^(i) becomes 
degenerate ast—^ ±oo. The implication is that the prob- 
ability for nonadiabatic transitions as a function of the 
pulse area is no longer expressed as a simple exponen- 
tial, Eq. ([3]), but rather by a sum of an exponential and 
an oscillatory term with an amplitude that vanishes only 
polynomially with the pulse area [111, [l^ • The exponen- 
tial term dominates for moderate pulse areas, whereas 
the oscillatory polynomial term dominates for large areas 
[lll[l3|- The border value of the area Ai,, where the expo- 
nential decline breaks down into slowly damped oscilla- 
tions, is proportional to the ratio R between the rms pulse 

area A 



' A"^ + A^ and the overlap area Ao (the area of 

overlap of the pump and Stokes pulses) Af, oc i? = A/Aq. 
The exponential decline of nonadiabatic transitions is fa- 
vorable for high-fidelity STIRAP because it allows one to 
achieve high fidelity with moderate pulse areas. This is 
turn implies a large value of the breakdown area Ah, so 
that the (slowly damped) oscillations emerge only when 
the infidelity is very low. Hence high-fidelity STIRAP is 
facilitated by asymmetric pulses, with longer outer tails, 
so that the ratio R (and hence the breakdown area Ai,) 
is large. These observations are further illustrated with 
the figures below. 

The above arguments suggest the following recipe for 
choosing Fit): 
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FIG. 1: Pulse shapes. 1: ideally optimized pulse shapes, 
which obey the optimization condition p7|l for all times t; 2: 
pulse shapes with a hypergaussian mask (|23|l . with n — 3, 
A = 4 and To = 2T; 3: Gaussian pulses (|27p . with pulse delay 
T = 1.2r. 



• Fit) should be a flat function during the time 
of overlap of the pump and Stokes pulses, during 
which the population transition takes place; 

• Fit) should have a sufficiently large width so that 
the rms-to-overlap ratio R is sufficiently large. 

In the examples below, as a mask Fit) we use the 
hypergaussian function. 



Fit) = e-(*/^«)' 



(23) 



where n = I corresponds to the Gaussian shape. For 
larger (positive integer) n the condition Fit) ~ const is 
fulfilled increasingly better in the overlap region. 

There is a some leeway in the choice of the function 
fit) as long as the adiabatic condition is fulfilled in the 
overlap region. We use 



fit) 



1 



-Xt/T " 



(24) 



We point out that other choices of the function /(t) are 
also possible. However, with the chosen method of opti- 
mization being based on the DDP approximation, which 
is valid only in the near-adiabatic regime, the function 
fit) has to fulfill the adiabaticity criterion 



(25) 



m\ « mi 



where i?(t) is the nonadiabatic coupling. Using Eq. ^5 
we obtain the following condition for the function /(t) 



\fit)\<t:no\Fit)\; 



(26) 
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FIG. 2: Deviation (infidelity) from complete population 
transfer vs the peak Rabi frequency for Gaussian (delay 
T = 1.2T) and DDP-optimized pulses (l23j (n = 3, A = 4, 
To = 2r and n = 1, A = 2, To = 2T). 



FIG. 3: Deviation (infidelity) from complete population 
transfer vs the peak Rabi frequency for DDP optimized pulses 
for different values of the width Tq of the hypergaussian 
(n = 3) mask F{t), Eq. ((23)) . 



hence f{t) should have a smooth time dependence in 
order to facilitate adiabaticity. Once in the adiabatic 
regime, however, the function f{t) does not affect the 
DDP optimization because it does not appear in the con- 
dition dni). 

Three different pulse shapes are shown in Fig. [TJ pulse 
shapes that obey the optimization condition p?|) for all 
times t, along with the more realistic example with a hy- 
pergaussian mask (|23|). that obey the optimization con- 
dition (|17p only in the region of overlap of the pump and 
Stokes pulses, and Gaussian pulses. 

np^noe-^'--/^y^\ r!, = r!oe-(*+^/2)/i^^, (27) 

where r is the pulse delay. We note that the pulse area A 
for the Gaussian pulses is almost identical with the DDP 
optimized pulses. In the following we demonstrate that 
the optimally shaped pulses are superior to the Gaus- 
sian pulses, even with optimized delay for the latter, in 
achieving a very high fidelity. 

B. Examples of ultrahigh-fidelity STIRAP 

In Figure ^ we plot the STIRAP infidehty, i.e. the 
deviation I — P3 from perfect transfer for the optimized 
pulses, described above, and compare these to the results 
for the traditional implementation of STIRAP with a pair 
of Gaussian pulses. The infidelity is shown as a function 
of the peak Rabi frequency. For Gaussian pulses, the 
pulse delay is chosen such that a nearly maximum fidelity 
is obtained. Despite the optimum delay, the Gaussian 
shapes do not allow us to reduce the infidelity below the 
limit 10~^ in the shown range of pulse areas (eventually, 
at very high pulse areas the infidelity drops below this 



limit). On the contrary, DDP-optimized pulse shapes 
(with n — 1, X = 2, Tq = 2T) break this limit even for 
small pulse areas. 

We note that the DDP pulses are not fully optimized 
due to the mask function F{t) vanishing at large times. 
Therefore the respective fidelity curve starts to oscil- 
late, which signals the occurence of nonadiabatic tran- 
sitions. However, the magnitude of these nonadiabatic 
transitions can be controlled by the width of the mask 
fmiction F{t): a larger width pushes these oscillations 
further down. Hence even for small pulse areas, the op- 
timally shaped pulses are far superior to the Gaussian 
pulses, as shown in Fig. [2l This tendency is also visible 
in Fig. [3] where the infidelity is plotted as a function of 
the pulse area for two different widths Tq of the mask 
function F{t), Eq. P5)) . By increasing the mask width 
Tq, the validity range of the adiabatic optimization con- 
dition (jl7p widens, and the optimized pulses approach 
the ideal DDP-optimized pulse shapes in Fig. [1] 

In Fig. |3]we compare the line profile as a function of the 
common detuning A from the intermediate level (see Eq. 
(121)) for the optimized and Gaussian pulse shapes. The 
time delay for the Gaussian pulses is numerically chosen 
for nearly maximum fidelity. It is known that a single- 
photon detuning does not affect the dark state (as long 
as two-photon resonance is maintained) 2]. Nonethe- 
less, adiabaticity deteriorates and the transfer efficiency 
decreases with increasing single-photon detuning. The 
robustness of the high fidehty STIRAP (i.e. STIRAP 
where the infidelity is below the limit lO^"*) against the 
single-photon detuning is much more pronounced for the 
optimized pulses. This feature is readily explained by the 
fact that the same pulse shapes optimize STIRAP both 
for A = and for large A, as discussed above. 
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FIG. 4: Line profile as a function of the common detuning A 
from the intermediate level for DDP-optimized (n = 3, A = 4, 
To = 2r) and Gaussian pulses (delay t = 1.2T) for a peak 
Rabi frequency f2 = 20. 



C. Relative error due to the RWA approximation 

The rotating wave approximation (RWA) is widely 
used whenever laser-induced excitations with optical fre- 
quencies to much larger than the Rabi frequency fl are 
considered. Typical optical (carrier) frequencies are 
Lu ~ 10^^ [s~^], while the typical Rabi frequencies are 



within the range f2 



10« 



10^[s~i] ta]. For two-level 



systems, perturbative inclusion of the counter-rotating 
terms results in the Bloch-Siegert shift of the eigenen- 
ergies by ^ fl'^/uj [13] ■ Within a typical pulse duration 
T - 10-6_io-8[s], this accumulates maximum deviation 
from the calculated probabilities of 



APe {T- 



IQ- 



10 



-16 



(28) 



It should be also noted that the results in RWA are exact 
in case of circular polarization Q. 



OPTIMIZATION OF HADAMARD GATES 
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FIG. 5: Raman line profile for the DDP-optimized (n = 3, 
A = 4 To = 2T) and Gaussian pulses (delay r = 1.2) for a 
peak Rabi frequency 57 = 20. 



Fractional STIRAP (f-STIRAP) is a variation of STI- 
RAP, which creates an arbitrary preselected coherent su- 
perposition of states ipi and tp3, 



"if = ipi cos a — -03 sin a. 



(29) 



As in STIRAP, the Stokes pulse precedes the pump pulse, 
but unlike STIRAP, where the Stokes pulse vanishes first, 
in f-STIRAP the two pulses vanish simultaneously, while 
maintaining a constant ratio of their amplitudes [13j , 



^ — oo< — t ^r>\i) t — >oo , 

< — „ , , — > tana. 



A convenient realization of f-STIRAP reads [ij 



(30) 



(31a) 



/2f/T^ 



-(t-T/2y/T^ 



cos Q;(»lb) 



The DDP-optimized pulses, in analogy with the full STI- 
RAP, read 



It is known that the transfer efficiency for STIRAP is 
much more sensitive to a detuning from Raman resonance 
than to the single-photon detuning . Figure [5] shows 
the Raman line profile for DDP-optimized and Gaussian 
pulses. As before, the time delay for the Gaussian pulses 
is numerically chosen such that a nearly maximum tran- 
sition probability P3 is obtained. The optimized pulse 
shapes are far superior to the Gaussian pulses and al- 
low one to maintain a high fidelity over a wide range of 
two-photon detunings. 



flpit) = noF{t)sin[af{t)], (32a) 
Qsit) - r!oF(i)cos[a/(i)], (32b) 

where f{t) is again an arbitrary function that satisfies 
condition (ffb)). It is easy to verify that these pulse shapes 
satisfy the f-STIRAP condition For half-STIRAP, 

when an equal coherent superposition of states ipi and ip^ 
is created, we should have a — n/A. This superposition 
corresponds to the Hadamard gate, which is one of the 
fundamental gates in quantum information processing. 

Figure [S] compares the transfer efficiency of f-STIRAP 
for DDP-optimized pulses and for pulses given by Eqs. 
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FIG. 6: Transfer efficiency for fractional STIRAP vs the peak 
Rabi frequency for DDP-optimized pulses (with n = 3, A = 4, 
To = 2T and n = 1, A = 2, To = 2r) and Gaussian pulses 
(pTjl (with r = 1.4T). 

([3T|) . The pulse delay for the pulses ((3T|) is chosen numer- 
ically for nearly maximum fidelity. In Fig. [6]the infidelity 
is defined as 

1 - I < inal I ^ desired > P, (33) 

where ^desired is given by Eq. (I29|) at a = 7r/4 and 
^ /inai is numerically calculated for both DDP-optimized 
and Gaussian pulses. It is important to note that despite 
the delay optimization, the pulses (|3T|) do not allow re- 
duction of the infidelity below a certain limit. Due to 
the mask function F{t), the DDP-optimized pulses led 
to an oscillatory behaviour of the infidelity. However, 
the oscillatory regime and hence the fidelity depend on 
the width To and can be controlled. In particular in the 



regime of small pulse areas, the DDP-optimized version 
of f-STIRAP (with n = 1, A = 2, To = 2T) is far superior 
to the traditional Gaussian pulses as shown in Fig. [6l 



VI. CONCLUSIONS 

We have proposed an optimization of the fidelity of 
the STIRAP technique, which uses the DDP approach 
to minimize nonadiabatic losses. The rationale for this 
is the reduction of STIRAP from three to two states ei- 
ther: on exact single-photon resonance or for large single- 
photon detuning. Interestingly, the optimized pulse 
shapes are the same in both regimes, which makes DDP- 
optimized STIRAP very robust against variations in the 
detuning. We have demonstrated with numerical sim- 
ulations that the fidelity of this optimized STIRAP can 
reach very high values, with an error well below the fault- 
tolerant QIP limit of lO"**, which is very difficult to reach 
by optimizing the pulse delay with the usual Gaussian 
pulse shapes (or other symmetric pulse shapes, such as 
hyperbolic-secant). The proposed optimization is of po- 
tential importance for QIP because it supplements the 
robustness of STIRAP against parameter variations with 
an ultrahigh fidelity of gate operations. We furthermore 
emphasize that a similar optimization method could also 
have a significant impact on the vacuum-stimulated Ra- 
man scattering in cavity QED [l^ . 
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